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Abstract 

We experimentally investigate the steady states of two granular assemblies differing in their material prop- 
erties and allowed to exchange volume with each other under external agitation in the vicinity of their jamming 
transition. We extract the statistics of various static and dynamic quantities, and uncover a materials-independent 
relationship between the average packing fraction and its fluctuations. This relationship defines an intensive 
parameter which decouples from the volume statistics, and remarkably takes the same value in both subsys- 
tems. We also observe that an effective diffusion coefficient also takes the same value in each subsystem, even 
as the structural relaxation time increases over several orders of magnitude. These observations provide strong 
constraints on the eventual estabUshment of a granular equation of state. 

1 Introduction 

Statistical mechanics has proven to be one of the most powerful tools for understanding transitions between states 
in physical systems. This predictive power has been successfully extended to gently perturbed systems through 
linear response theory. ' Such an approach emphasizes the significance of intensive state variables, analogous to 
temperature or pressure in equilibrium systems, but potentially extends their relevance to the wider spectrum of 
out-of-equilibrium systems. However, this formalism is not valid for systems driven far from equilibrium, even 
when in well defined steady states. It remains an open question whether or not state variables can be defined in a 
generic way for non-equilibrium steady states. The hope is that such state variables would carry predictive power 
regarding the system's response to a change in the driving parameters, or its behavior when put in contact with 
another system. 

A key situation in which such tools are expected to be relevant is that of granular materials: systems composed 
of individual particles large enough to be insensitive to thermal fluctuations. Such materials are ubiquitous in 
natural and industrial contexts, and models of their dynamics would benefit from an improved understanding of 
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how to prepare and manipulate them in well-defined macrostates determined by a small set of control parameters 
analogous to thermal intensivities. Promising measures of the temperature-hke Edwards compactivity^'^ have 
been made in static granular systems^, as well as Kubo-type fluctuation/response pairs ^ leading to effective 
temperature measurements in dynamic granular systems ^"^ ^ . However, it remains an open question to what extent 
these temperature-like variables behave as true temperatures. 

For granular materials under strong agitation, the interactions of individual particles are dominated by binary 
collisions. Such dynamics are well-described by a thermal-like kinetic theory in which the velocity fluctuations 
play the role of temperature However, the denser states relevant to many important granular flows are dom- 
inated by multi-body interactions and such states remain difficult to model through kinetic theory approaches^"*. 
The phenomenology of these dense states is strikingly reminiscent of supercooled liquids 1^-20^ where classical 
statistical mechanics also loses validity due to a loss of ergodicity. Both the thermodynamic glass transition 
and the jamming transition in dense granular materials lack a fundamental understanding, and both have been 
hampered by the difficulty of making a relevant quantitative description of the state of the system. This universal 
phenomenology, embodied by a dramatic slowing down of the dynamics'^^ caging and memory effects 
is now referred to as the jamming transition, and is at the heart of a very active field of research which is the subject 
of this special issue. 

Thermodynamics is founded on, among other things, a meaningful equilibration of intensive state variables be- 
tween subsystems. Thus, if a thermodynamic-like description of granular materials is to be found, an informative 
first step would be to search for an analog of the zeroth law of thermodynamics. By putting two athermal systems 
in contact and allowing them to reach a steady state, one can seek relevant intensive parameters by discarding 
properties that fail to equihbrate. We experimentally conduct such an investigation, putting two model granular 
subsystems in contact. 

In our experiments, we aUow two adjacent monolayers of particles, differing in material properties, to ex- 
change volume with each other under external agitation and an overall constant- volume constraint. We study the 
steady-state properties of the two subsystems as a function of the overall average packing fraction in the vicinity 
of their jamming transition. Importantly, we uncover a materials-independent relationship between the average 
packing fraction and its fluctuations. This relationship defines an intensive parameter that decouples from the 
volume statistics of the subsystems, and remarkably takes the same value in both subsystems. This observation 
provides a strong constraint on the eventual establishment of a granular equation of state. Moreover, among 
various dynamical quantities studied, an effective diffusion coefficient is shown to take the same value in the 
"equilibrated" subsystems, while the structural relaxation times exhibit different behaviors. While a thermody- 
namic system would lead to temperature equilibration, the present situation leads to a decouphng of the mobility 
and the structural relaxation time, and only the former seems to yield relevant information about the properties 
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Figure 1: Apparatus schematic (1 m x 2 m), with particles drawn to scale. Filled particles are taken from an 
experimental configuration and indicate the size of the imaged region. A between and 1 measures the fractional 
position of the piston. 

shared by subsystems in steady states. Remarkably, this decoupling of the dynamics occurs in the vicinity of the 
static random loose packing. 

2 Experimental setup and protocol 

We study the equilibration of two adjacent dense bidisperse layers of disks on an air table, as sketched in Fig. 1. A 
piston separates the two subsystems, and is constrained to move along the long axis of the air table by protruding 
feet which impede rotation. On each side of the table, the system is prepared with an identical number of particles 
A^, with the ratio of small to large particles fixed at Ns = 2Nl- This ratio suppresses crystallization and provides 
roughly the same volume of large (diameter = 86 mm) and small particles (ds = 58 mm). The particles within 
each subsystem have a different set of material properties: on Side 1, the particles have restitution coefficient 
£1 = 0.5 1 ± 0.07 and friction coefficient jUi = 0.85; on Side 2, the particles have £2 = 0.33 ± 0.03 and jU2 = 0.5. 
Restitution coefficients were measured from isolated binary collisions; friction coefficients are nominal values 
from the literature. Both sides utilize standard plastic Petri dishes as the particles, with the difference in particle 
properties achieved by encircling the particles on Side 1 with a rubber band. Because the sides of the particles 
slope inwards, the thickness of the rubber band does not significantly change the radius of the particles; the mass 
of the particles on Side 1 is increased by 7%. 

The aggregate rearranges via an array of sixty electromagnetic bumpers which form the walls of the system. 
These bumpers are triggered pairwise: bumpers facing each other in the system fire at the same time in order 
to prevent net momentum and torque injection. These pairs are triggered randomly via a pre-generated random 
sequence. Four pairs of bumpers are randomly fired every 0. 1 second, and travel 1 cm into the granular pack; the 
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total time during which the bumpers stay in their forward position is approximately 0.1 sec. 

To quantify the long-time mobiUty of the particles, we take images at a frequency which is low compared 
to the energy injection timescale, making usual tracking techniques inoperative. Therefore, we have developed 
a tracking method which identifies each particle by a unique tag. Each particle is marked with a 3 x 3 array of 
colored dots which encodes two copies of a 4-bit, 4-digit identifier, plus an error-correcting bit. The particles are 
located by their circular rims and their identities are established using the tags, allowing their positions adjacent 
image frames to be connected into trajectories. We monitor the positions of the piston and the inner 75% of the 
particles with a CCD camera mounted above the apparatus; we obtain a minimum of 10"* configurations for each 
experiment. 

To imderstand the equilibration of the system on its approach to janoming, we perform experiments at increas- 
ing values of N while holding all other variables constant. We adjust A'^ in increments of 2 small and 1 large 
particle on each side, so that the two sides always have the same A'^. The average packing fraction ^ calculated for 
the entire system is given by 



as well as by the ratio of the total area of A'^ disks to the total area of the air table. As we adjust N from 183 to 204, 
this corresponds to ^ = 0.768 to 0.818. 

3 Macroscopic observables 

First we focus on the macroscopic state of the system, namely the volume occupied by each of the subsystems. 
Since the overall volume of the system is conserved, we only need to consider the position of the piston scaled by 
the length of the cell < X < 1, where the origin of the axis is taken on the side occupied by the particles that 
have been circled by a rubber band, referred to as Side 1 in the following discussion. 

We first consider the time signal A (?) for a given value of A^, an example of which is displayed in Fig. 2a. After 
a transient time Tr, the signal becomes stationary and the system has reached steady state. This state is observed 
for all reported values of A'^. Note that Tr is comparable to the a-relaxation time of the assembly, i.e. the time 
it takes for the particles to diffuse further that their size^^. We hence consider the probability distribution of A in 
this stationary regime and we find Gaussian statistics even for the largest A'^ (densest ^); sample distributions are 
shown in Fig. 2b. 

These statistics are thus characterized by the first two moments. We observe that the average piston position 
(X) = {X {t))t is observed to be < ^ for all values of A^: Side 1, containing the particles with rubber bands, 
occupies less volume than the side without. We conclude that this systematic deviation from equal volumes 
originates from the difference in material properties of the particles. Furthermore, and rather surprisingly, (A) is 
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Figure 2: (a) Scaled piston position A (?) for = 201 {(j) = 0.8 12). (b) Sample probability density functions of X 
(and Gaussian fits) for A? = 183, 192 and 201 = 0.766, 0.789, 0.812). 
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Figure 3: (a) (A) (Side 1) as a function of A^, with bars representing the standard deviation. For some values of A^, 
the experiment has been repeated several times starting from different initial configurations, in order to evaluate 
the errors; these measurements are shown as smaller symbols, (b) Standard deviation 5X as a function of A^. 
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not a monotonic function of A'^: it increases until it reaches j for AT* = 195 and then decreases for larger values 
of A'^. At this crossover, the packing fraction is equal on both sides, <j)* = 0.798. The behavior of the system is 
qualitatively different on each side of this transition point. For values of N <N*, the particles interact mostly 
through short binary colUsions and the system looks liquid-Uke. For N>N*,the dynamics is dominated by multi- 
body interactions and becomes much slower. In watching the system, intermittent force chains'" are be observed 
to carry the motion imposed by the boundaries into the bulk of the system, and some particles stay in contact 
with one another on time scales larger than the typical collision time. However, no permanent stress seems to be 
trapped in the system until it jams: this occurs at the highest explored value of A' as wiU be quantified below. 

Interestingly, this transition at ^* = 0.798 is close to, but perhaps slightly below, independent measurements 
of random loose packing for these particles: <j)^ = 0.807 ± 0.010 and (j)^"' = 0.812 ± 0.006. To determine (j)'^^ 
for each of the two types of particles, we placed the table at a 0.3° angle, and rained down single particles from the 
high end to the low end to create the loosest packing accessible to us. These values are upper bounds for random 
loose packing, for our protocol does not guarantee that we reach the loosest possible stable packs. 

The standard deviation of these distributions, 5A = (A^) — (A)^, exhibits a simpler behavior: it is a decreas- 
ing function of A', as can be seen in Fig. 3b. This is to be expected from the fact that as A' increases a smaller 
amount of free volume is available for the motion of the piston. Notably, this quantity does not exhibit any special 
feature in the vicinity of A'* . Finally, the magnitude of dX at the highest explored value of A^ closely matches the 
volume fluctuations imposed by the bimipers. This mean that this macroscopic quantity becomes trivial when the 
system approaches its jammed state, merely reproducing the microscopic fluctuations imposed by the injection 
mechanism. However, even in the densest state explored the system is not strictly rigid: particles continue to 
rearrange due to the impacts from the boundaries. 

4 Microscopic observables 
4.1 Statics 

We now turn to the microscopic information provided by particle-tracking. In order to extract a well-defined, size- 
independent microscopic quantity from the packing fraction field, we analyze the spatial and temporal fluctuations 
of the average local packing fraction <j) over windows of increasing size, where <j) is defined as the ratio of the 
number of pixels occupied by the particles to the total number of pixels in a given region. For each subsystem, we 

measure (j) over boxes of size L ranging from a few = ^d^ up to half the system size. For L > Ir^, (j) {!-■) converges 
to a constant. Moreover, the variance (5^^) scales approximately as L^^, as shown in Fig. 4. Such scaling behavior 
is expected from the central limit theorem provided there are no long-ranged spatial correlations in the packing 
fraction field. Since the packing fraction does not exhibit such correlations, it is suitable for a thermodynamic-like 
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analysis despite the relatively small number of particles. Furthermore, we obtain an L-independent measure of the 
variance of (j) by averaging (50^)o = {5<j)'^)L? over all L > 2rs. 

These two packing fraction statistics are plotted as a function of N in Fig. 5: {(j)) is higher for Side 1, in 
agreement with our macroscopic measurement. Similarly, at N = 195 the packing fractions are equal {{(j)) = 
^ = 0.798 for both sides). Different values of the normalized ^ fluctuations (50^)o are observed on the two 
equilibrated sides, indicating that the statistics of this quantity depend on the material properties. Moreover, the 
fluctuations on each side are found to be decreasing functions of A^, as shown in Fig. 5b. As for the fluctuations of 
A, this can be understood by the fact that as N increases, the amount of free volume to be distributed among the 
particles decreases. Finally, the discrepancy in the packing fraction statistics between the two sides gets larger on 
approach to jamming at large ^ (large AO, especiaUy as the number of particles is increased above A'^*. 

To examine the relevance of these two state variables, we plot the dependence of (50^)o on (0), parametrized 
by A'^ (see Fig. 5c). Remarkably, the data from both sides fall onto a single master curve: the fluctuations in local 
are insensitive to the material properties of the particles. Within the explored range of A'^, this master curve is 
approximately linear and the extrapolation would intersect (50^)o = at (0) = 0/ = 0.845, which is compatible 
with the rigidity transition reported in other bidimensional systems ■ ' . This result therefore stands as a starting 
point for a granular equation of state. It is notable that, as in Fig. 5b, the transition at A'* = 195 does not appear as 
a feature. 

4.2 Dynamics 

We now turn to the microscopic characterization of the dynamics in our system as the system transitions from 
a liquid-like to a solid-like state. To monitor this transition, we use particle-trajectories to compute the average 
diffusion distance as a function of lag-time T, defined as 



where r, [t) is the vector position of particle i at time /. Averages are computed over the ensemble of all particle 
trajectories at all times. 

Fig. 6 shows how the r.m.s. displacement g{x) and the corresponding particle trajectories vary as increasing 
A'^ brings the system towards jamming. AtN= 183 (low ^) a particle explores a region several wide, and a 
saturates at long time scales due to finite system size. On the other hand, atN = 201 (high ^), a exhibits sub- 
diffusive behavior at short T, and caging effects are significant. As the dynamics of the system slow down at large 
N, we increase the duration of the experiment from twenty hours up to fifty-five hours for the densest packing in 
order to sample a significant set of configurations. However, even at this extended acquisition time, the diffusion 
length atN = 201 barely passes Ir^. From the glass transition point of view, the system is completely jammed. 
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Figure 4: Variance (50^) of the packing fraction measured within squares of size L for three values of A^; dashed 
Hnes are {5^^)o/L^. 
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Figure 5: (a) Packing fraction {(j)) and (b) normalized packing fraction fluctuations (50^)o as a function of number 
of particles on each side, (c) Data from parts (a) and (b) combined on a single plot without regard for A^. Ellipse 
encloses two points obtained from Side 1 (red) and Side 2 (black) for a single run with = 186. Dashed line is a 
linear fit to data from both sides, showing an intercept at = 0.848. 
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Figure 6: (a) Trajectories of a single particle (from Side 1) for a duration of 10 sec. (b) Diffusion length a as a 
function of lag-time T. Red lines are for Side 1; black lines are for Side 2. Dashed line is a °^ to guide the eye. 
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Figure 7: Local exponent v as a function of lag-time T for = 183, 192 and 201. The thick region indicates the 
averaged maximum used to obtain a diffusion constant D from Eq. 3. 

In order to extract meaningful information from (t(t), including establishing whether a truly diffusive regime 
is ever reached, we compute the local slope 
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d{logx) ' 

over a range of T that allows for a reasonable estimation of the local slope. Interestingly, for all studied values of 
this local exponent exhibits a maximum, as illustrated in Fig. 7. This maximum is consistently close to j, as 
can be seen in Fig. 8b. However, the exponent V is larger for Side 1 (rubber bands) than Side 2 (bare) for < A^*. 
Above A^*, Side 1 changes from super-diffusive (v > j) to sub-diffusive (v < j) behavior. In contrast. Side 2 is 
sub-diffusive for all values of A^. 

We have extracted the characteristic time and length scales Tm and Om at which this maximum is reached, 
plotted in Fig. 8. For both sides of the system, the time scale at which maximally-diffusive behavior is reached 
rises by two orders of magnitude as A^ approaches the jammed state, and the length scale falls. These quantities 
indicate the occurrence of a dynamical crossover at A^* since Gm, Tm, and v all display different behavior above 
and below this value. Side 1 is more mobile for N < N* and Side 2 is more mobile above. Interestingly, this 
transition occurs when Gm is close to the mean diameter of the grains {Gm = 2ri). This suggests that the local 
rearrangement mechanism governing the dynamics at the particle scale changes near A^*. 

We use the behavior near {Gm,'^m) to extract quantities which quantify the diffusive behavior of the system. 
The structural relaxation time Ta is obtained from the time at which C7(Ta) crosses r^, meaning that the particle 
has diffused beyond its radius. A diffusion coefficient D is estimated by solving 
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for D. Since the dynamics of the system are not properly diffusive, this diffusion coefficient should be interpreted 
merely as a measure of the displacement amplitude at the time of maximum mobility. Astonishingly, as can be 
seen in Fig. 8, this definition for D provides a dynamical parameter that takes approximately the same value in 
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Figure 8: Maximum exponent Vm, and corresponding values of Gm and Tm as a function N for the two grain 
species, (c) Diffusion coefficient D and (d) structural relaxation time Ta as a function of A^. 



each subsystem regardless of A^. 

Both D and Ta indicate a slowing down for increasing A^, as shown in Fig. 8. On both sides, D decreases mono- 
tonically to nearly zero for our last data point, indicative of the system's cessation of rearrangement (jamming). 
Meanwhile, as is usually observed in the vicinity of the glass/jamming transition, the structural relaxation time for 
both species of particles soars over multiple orders of magnitude. The plot of Ta (N) shows that both above and 
below A^*, Side 1 (rubber band particles) exhibit longer relaxation times than Side 2. As for v, Gm, and Tm, the 
A^* crossover point is special in that these quantities each take the same value on both sides. No dramatic changes 
take place for D{N) near A^* : it is remarkably equal in each subsystem for all values of A^. 

Since the timescale Ta is observed to be sensitive to A^*, we investigate the extent to which the packing fraction 
is relevant to the dynamics. Therefore, we plot both D and Ta as a function of {(p) in Fig. 9. As in Fig. 5, pairs of 
data points from the same run are plotted at their corresponding values of (0). This parametrization predicts that 
D{^) vanishes at different values for the two different particle types, with ^j'-* < In examining Ta(0) for the 
two sides, there is a common (material-independent) branch of the curve for (0) < (j)* which increases gradually. 
The value of {(p) at which Ta leaves this branch differs for the two types of particles, with Side 1 again having its 
transition at lower {(j)). This transition is associated with local rearrangement mechanisms, and likely corresponds 
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Figure 9: (a) Diffusion coefficient D and (b) structural relaxation time Ta as a function of {(p) (without regard for 



to the packing fraction at which single-particle rearrangements are no longer possible 



5 Discussion and Conclusions 

We have explored the stationary states of subsystems able to exchange volume, following them from a liquid-like 
to a jammed state. In order to examine the extent to which state variables and the density of states are useful tools 
to understand this transition, we would need to understand the role of energy, entropy, and our ability to identify 
the corresponding microstates. An analysis based on the generation of entropy through the braiding of trajecto- 
ries ^^'^^ establishes the chaotic nature of the internal dynamics in our system. This, together with the stationarity, 
ensures the existence of an invariant measure, namely the SRB measure from which our temporal averages 
are guaranteed to correspond to ensemble averages. This provides a strong basis for the statistical analysis pro- 
vided above and lends a thermodynamic-like status to the quantities extracted from our equilibration experiments. 
Within such a framework, our results suggest further interpretation. 

In Fig. 5, we have uncovered an unexpected one-to-one relationship between the average local packing fraction 
(0) and its fluctuations (50^)o. This striking relationship, in which both and the material properties collapse 
on a single curve, should be a prediction from an as-yet unknown equation of state. Interestingly, the function 
does not appear to be sensitive to the transition at 0* « (j)'^^^ which is a prominent feature in the equilibration 
measurements of the piston position (see Fig. 2). 

To make use of the relationship, it is tempting to define a variable 



(5) 



which is intensive, ^-independent, and takes the same value in both subsystems for all A'. However, the value of 
Q, which we conjecture to be related to the energy injection and/or dissipation rates in the system, is not sufficient 
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to specify the equilibrated volumes. This suggests that another parameter is needed to fuUy characterize the state 
of a dynamically evolving dense granular pack, much as both temperature and pressure equilibration are needed 
to solve the equivalent problem within classical thermodynamics. The present experiment does not allow us to 
measure the instantaneous kinetic and potential energies, but future studies of particle kinetics would allow further 
investigation of such parameters. 

Diffusion measurements provide an alternative framework in which to understand the behavior of our system 
and hence the equilibration process. We have observed that the piston position X, relaxation time Ta, and the 
maximum mobility length and timescales (cTm, Im) are sensitive to the transition at (^* (see Fig. 2) which is close 
to the static 0^^. These changes in the dynamics suggest that the transition corresponds to two quahtatively 
different dynamical regimes. For ^ < 0*, the particles primarily interact through short binary collisions and 
there is weak coupling between the particles. Because the system is colUsion-dominated, it is likely that the 
restitution coefficients of the particles are an important consideration in the dynamics. If so, it would explain the 
observation that for Side 1 (larger coefficient of restitution), A increases with N. For ^ > 0*, long multibody 
interactions dominate the dynamics, and we visually observe signs of intermittent force chains^" which strongly 
couple both particle-particle and particle-bumper interactions. It is likely that the higher friction coefficient on 
Side 1 significantly increases the structural relaxation time of the corresponding assembly, and the dynamics are 
more sensitive to the frictional properties of the particles. Interestingly, all measured quantities are approximately 
equal in each subsystem at N*. 

As the system approaches ^j, the persistence time of the contact network increases, and the internal response 
of the particles is more correlated. For the largest studied value of N, the packing fraction and overall volume 
fluctuations are are imposed by the boimdaries. Hardly any rearrangement occurs on experimentaUy-accessible 
time scales, but extended internal vibrations can be observed: the behavior of the system is thus comparable to 
that of a solid. 

From a thermodynamic point of view, it is remarkable that the effective diffusion coefficients defined above 
are the only measured quantities to take the same value in the equilibrated subsystems, whereas the relaxation 
times Xa become very different from each other beyond A^*. As a matter of fact, in equilibrium systems, the 
Stokes-Einstein relation states that the product of the diffusion coefficient and the viscosity (or, equivalently, the 
relaxation time), is proportional to temperature ^ Hence, only the product Dt^, and not D itself, is expected to 
equilibrate. Here, we find that there is a decouphng between the relaxation time and the diffusion coefficient 
beyond N*, which is reminiscent of the decoupling reported in the vicinity of the glass transition in supercooled 
Uquids^^'^*, the latter equiUbrating while the former emphasizes the difference in material properties between the 
two sides. This phenomenon is yet another common feature between dense granular materials and supercooled 
Uquids, and seems to indicate that fluctuation/response based temperatures might be relevant below (^^^. 

Thus, we observe that both packing fraction (or free volume) and material properties play a key role in de- 
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termming the state of the system. In static granular systems, the Edwards ensemble^ provides a framework for 
imderstanding how volume plays a central role in describing the density of states and associated Boltzmann-like 
thermodynamic properties. While our system is dynamic, it is interesting to ask to what extent the equilibration 
of volume between two subsystems exhibits features reminiscent of the Edwards formaUsm. The overall vol- 
ume of our system is a conserved quantity, the dynamics explore a set of configurations (although not strictly 
mechanically-stable), and we achieve stationary distributions of system properties such as {(p). Were volume the 
dominant state variable setting the density of states for this system, we would have observed A = j, independent 
of A^. Instead, the material properties of the particles play an important role in dynamically determining the state 
of the system, and we observe that in general A ^ j. 

In this paper, we have experimentally studied the steady state statistics of two model granular subsystems 
differing in their material properties put in contact through a mobile wall. We have uncovered several features 
that constrain any specification of a granular equation of state. First, most static and dynamical quantities are 
affected by a dynamical crossover which occurs at a packing fraction close to the (static) random loose packing. 
We interpret this crossover as separating the weak and strong coupUng regimes in the behavior of the system. 
However, both the macroscopic volume fluctuations and the diffusion coefficients are insensitive to this transition. 
Second, we have found that the packing fraction fluctuations are in a one-to-one correspondence with the average 
packing fraction, independent of the material properties of the grains. These fluctuations drop to zero on the 
approach to jamming. Finally, the diffusion properties of the particles reveal a glassy feature known as decoupling 
which arises above random loose packing. However, in spite of the fact that the structural relaxation timescale 
becomes both large and material-dependent for dense systems, the effective diffusion coefficient nonetheless takes 
the same value in each subsystem. This suggests that this coefficient is an important intensive parameter describing 
the state of the system. A theoretical understanding of these last two findings could come from the shrinking 
density of states on the approach to jamming. It is intriguing that the packing fraction fluctuations (geometry) 
suggest material-independence while the relaxation timescales (dynamics) suggest material-dependence. 
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